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Furnishing the Galaxy with Pulsars 
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Australia 

Abstract. The majority of pulsar population synthesis studies performed to date have focused on isolated pulsar evolution. 
Those that have incorporated pulsar evolution within binary systems have tended to either treat binary evolution poorly of 
evolve the pulsar population in an ad-hoc manner. Here we present the first model of the Galactic field pulsar population that 
includes a comprehensive treatment of both binary and pulsar evolution. Synthetic observational surveys mimicking a variety 
of radio telescopes are then performed on this population. As such, a complete and direct comparison of model data with 
observations of the pulsar population within the Galactic disk is now possible. The tool used for completing this work is a 
code comprised of three components: stellar/binary evolution. Galactic kinematics and survey selection effects. Here we give 
a brief overview of the method and assumptions involved with each component. Some preliminary results are also presented 
as well as plans for future applications of the code. 
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POPULATION SYNTHESIS 

We evolve n systems by first selecting the initial param- 
eters of each single star (mass) or binary (masses, pe- 
riod) from appropriate distributions, e.g.- stellar masses 
drawn from the Kroupa, Tout & Gilmore (1993) initial 
mass function. Limits must also be chosen where we typ- 
ically take 5 - SOMg, and 0.1 - 4-QMq, for the binary pri- 
mary and secondary masses. The algorithm loops over 
the n systems in binpop which follows the stellar/binary 
evolution. Any systems of interest that are produced (in 
this case pulsars) are passed to binkin which follows the 
kinematic evolution of each system. The resulting popu- 
lation is then 'surveyed' in binsfx which produces syn- 
thetic data for analysis. 



BINPOP: BINARY AND PULSAR 
EVOLUTION 

This component evolves single and binary stars from the 
zero age main sequence through to some given time, t^ax- 
It uses the binary evolution code of Hurley, Tout & Pols 
(2002) with the addition of pulsar evolution. Some fea- 
tures included in the binary evolutionary algorithm are 
tides, mass transfer (examples are stellar winds, roche- 
lobe overflow and common envelope), magnetic brak- 
ing, orbital gravitational radiation, mergers and super- 
nova velocity kicks. Features particular to neutron star 
(NS) include, standard pulsar magnetic field decay, ac- 
cretion induced field decay and spin-up, magnetic brak- 
ing, propeller effects and death lines. In this body of 
work we focus on a particular model which makes the 



following assumptions and produced flie pulsars shown 
in Figure 1. The model evolves n = 10^ binary sys- 
tems, where each system is evolved to a f„„„ drawn at 
random from between the limits of and 10 Gyr If a 
NS forms it's initial spin period Pq and magnetic field 
Bo are given by Pq = 0.02 + 0.156{Vkick/Va) [s] and 
Bo = 5x1011 +4xl0i2 X (Vkick/Va) [G], where VkM is 
randomly drawn from a Maxwellian with standard devia- 
tion Va = 190 km/s. In terms of NS magnetic field decay 
we assumed the standard decay time constant , Tg = 1000 
Myr while the accretion induced decay scaling param- 
eter is, k = 7000. We then ignore pulsars whose mag- 
netic field decays to the limiting value of 5 x 10^ G. If 
a NS is accreting we assume it is not visible as a pulsar 
and it is ignored, this includes those NS who were ac- 
creting material via Roche-lobe overflow within the pre- 
vious 1000 years. Finally, this model also includes the 
Harding, Muslimov & Zhang (2002) curvature radiation 
death line, beyond which pulsars are assumed not to emit 
beamed emission of sufficient strength to be observed. 



BINKIN: GALACTIC KINEMATICS 

In this component single stars and binary stars are 
evolved dynamically within the Galactic gravitational 
potential, which we model as given by Kuijken & 
Gilmore (1989). This is our first and simplest model of 
the Galaxy and includes the disk, nucleus and bulge. 
Binkin also takes into account supernovae, which may 
impart some momentum onto a single star or binary sys- 
tem or, in fact, disrupt a binary system. Every system 
has it's orientation randomised with respect to the galac- 
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FIGURE 1. Spin period against spin period derivative for 
each pulsar produced in the model described to the left. Each 
point is one pulsar, circles represent binary pulsars while stars 
are isolated pulsars. 



tic plane. The code has been tested and is reversible. 
We note that this is not an n-body code, each system is 
evolved independently. Each of the pulsars in Figure 2 & 
3 are evolved in binkin which moves them in the Galaxy 
to produce an 3-dimensional spatial distribution. In this 
model the stars are assumed to be born in a disk with a 
radius of 15 kpc and a height above the plane of 50pc. 



BINSFX: SYNTHETIC SURVEYS 

We now consider the final component of our work, 
binsfx, which uses both the pulsar data from binpop and 
the Galactic positional data from binkin to produce syn- 
thetic surveys. These synthetic surveys are based on the 
Parkes Multibeam (PMB), Swinburne intermediate lati- 
tude (SIL) and Swinburne extended (SE) surveys along 
with the Burgay et al. (Betal) survey. We also make use 
of the GMRT, which has not yet been used for an ex- 
tended and large pulsar survey project. In modeling the 
surveys it is important that we also take into account the 
associated pulsar selection effects. Within the model pre- 
sented here (Figures 2 & 3) we have considered beam- 
ing of the pulses away from our line of sight, scatter- 
ing of the image (multi-path propagation), dispersion (a 
lag between different wavelengths) and sampling of the 
pulse with a finite time constant. We also need to assume 
some form of the pulsar luminosity and we do this using, 

L = 0.1 X 10"-i^ x-P" X ^^ where a = -2.0 + O.lx, 
j3 = O.O + O.lj and j is a Gaussian random number. 
We assume a spectral index of the form — O.Sj. Once 



we check whether each pulsar produced above (depicted 
by Figure 1) is within the correct region of the Galaxy 
for the particular survey and is close enough that our 
assumed telescope sensitivity is less than the calculated 
pulsar luminosity (at the solar system) then we may start 
to produce plots like Figures 2 & 3. Note that Figure 2 
is a reproduction of Figure 1 after passing the population 
through binkin and binsfx. We may now compare these 
to the observations. For our simulated surveys of PMB, 
SIL, SE and Betal we find 1 180, 105, 54 and 99 pulsars 
respectively. With the above assumptions we predict that 
the GMRT could see 366 pulsars. 



FUTURE WORK 

The results shown here are a guide to what we can now 
achieve when modelling pulsars within the Galaxy. Be- 
ing able to directly compare the results with observa- 
tion, although not new, is a step forward when simulat- 
ing with such detailed stellar and binary evolution codes. 
For example, the number of pulsars the above simu- 
lated surveys 'see' are obviously only estimates, how- 
ever, through comparison with pulsar surveys we are now 
able to better constrain the complete underlying pulsar 
luminosity law. This work will also help in constraining 
binary and pulsar theory by allowing us to compare the 
spatial distributions of differing system types (for exam- 
ple see Figure 3) and the numbers and birth rates of these 
systems as directly compared to observations. Note that 
the binary evolution algorithm is not restricted to pul- 
sars and any star/binary of interest can be investigated. 
This work will also help predict what next generation 
telescopes, such as the Square Kilometer Array, should 
see. 
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FIGURE 2. Synthetically observed pulsars from the Galactic pulsar population shown Figure 1. Each point is one pulsar, crosses 
are isolated millisecond pulsars while circles are binary millisecond pulsars. 











- 


















. 


















^ 


,^ 


^ '■* 






' 










r 


|fe 


s 














1 


mm-.-'r' 


■" 














»' a ; •> * 














" ?• .'■••',•■7. 




*^ X ". 








o 






. ' ■•' • 




a 




' 








X 













X [t:pc] 



FIGURE 3. The Galactic x-y positions of those synthetically observed pulsars shown in Figure 2. Each point is one pulsar, 
crosses are isolated millisecond pulsars while circles are binary millisecond pulsars. 



